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Single-molecule spectromicroscopy: the door into sub-diffraction refractometry 
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We suggest a novel approach for probing of local fluctuations of the refractive index n in solids 
by means of single-molecule (SM) spectroscopy. It is based on the dependence Ti(n) of the effective 
radiative lifetime Ti of dye centres in solids on n due to the local-field effects. Detection of SM zero- 
phonon lines at ultra-low temperatures gives the values of SM natural spectral linewidth (which is 
inverse proportional to Ti) and makes it possible to reveal the distribution of the local n values 
in solids. Here we demonstrate this possibility on the example of amorphous polyethylene and 
poly crystalline naphthalene doped with terrylene. 
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One of the most important characteristics of a mate¬ 
rial that determine many of its macroscopic properties is 
the refractive index n. It is included as a parameter in 
classical equations (e.g., Fresnel’s and Maxwell’s), linked 
to dielectric permitivity and magnetic permeability, and 
known to depend on temperature, pressure and wave¬ 
length [I[. In recent years a special attention has been 
attracted to the refractive index in relation to metama¬ 
terials, which are assumed to have the negative n [i,i- 

Since Ernst Abbe invented the refractometer in 1874 
0 the most precise techniques for measuring n have 
been based on the principle of the total internal reflec¬ 
tion (TIR) on the interface of two media. There are also 
some other practical techniques for direct and indirect 
measurements of n (goniometry, elli^ometry, interfer¬ 
ometry, frustrated TIR, see, e.g. however, most 

of them require averaging over a large macroscopic vol¬ 
ume of the sample. As a result the question of the lo¬ 
cal fluctuations of n in real solids remains unanswered. 
Moreover, in view of the growing interest to nanotechnol¬ 
ogy and nanophotonics, it is important to know whether 
there is a relation between n and micro- and nanoscopic 
structure of solids. In order to give the answer some 
special instruments for probing the local n values on the 
sub-diffraction (nanometer) level must be developed. 

The topical papers describe a great number of very 
promising “nanoinstruments” for probing of various lo¬ 
cal parameters of a sample based on single quantum light 
emitters embedded into the sample’s material. For exam¬ 
ple, single-molecule (SM), single quantum dot (QD) and 
single nanocrystals spectroscopy and imaging provides 
an outstanding basis for various kind of nanoinstrumen¬ 
tation, like nanoscale thermometry nanodetectors of 
acoustic strain 0, single charges detectors [ll|, instru¬ 
ments for probing of low-energy excitations in disordered 
solids |I2l4l7j. etc. Besides, we see a possibility to use 
SM spectroscopy for probing the local-field effects in dye- 
doped solids. It is well know that the local-field effects 
cause the dependence of the exited state lifetime of a light 
emitter on the refractive index of the hosting medium 
Ti(n) Thus, the study of photo-physical proper¬ 


ties of a single dye center is the way to probe the local n 
values. For example, in [ 2 ^ such an idea was realised by 
analysing the fluorescence decay of QDs. 

In this work we propose a method for the probing of the 
local (nm-scale) fluctuations of n in solids by the analysis 
of zero-phonon spectral lines (ZPL) of single impurity 
dye-molecules at ultra-low temperatures (Fig. [T]). 


ZPL, which corresponds to aj)urely electronic transi¬ 
tion in an impurity molecule [271, is a unique source of 
information about dye-matrix interactions [l2|, [l3|. ZPL 
parameters (frequency, intensity, width, etc.) are very 
sensitive to the local environment of the corresponding 
chromophore SM. This fact makes SMs good candidates 
for spectral probes to obtain data on the structure and 
the internal dynamics of solids. 


As we know from numerous studies, the temperature 
T dependence of the homogeneous spectral width of ZPL 
is determined by three main contributions 12 fH, i.e.: 


^ZPl{T) = To + /SM e-tunn{T ,tm) + /ST e-phon{T), (I) 



FIG. I. The illustrative sketch of the developed approach for 
the probing of local values of the refractive index n in solids. 
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where the natural (lifetime limited) linewidth Fq is 


ro = 


1 

27rTi’ 


( 2 ) 


the ZPL broadening due to the interaction of electronic 
transitions in impurity molecules with tunneling excita¬ 
tions in a matrix is 


~ r“ln(tm); (3) 


the ZPL broadening due to the quadratic electron- 
phonon interaction, which in the simplest case of inter¬ 
action with a single quasilocalized low-frequency vibra¬ 
tional mode (LFM) is expressed as 


e—phon{^) — 


ex.p{—AE/kT) 

[1 — ex.-p{—AE/kT)]‘^ ’ 


( 4 ) 


In these equations tm is the total time of measurement, 
1 < a < 2 and its exact value depends on the dye-matrix 
system, w is the constant of quadratic electron-phonon 
coupling, AE is the energy of LFM. 

ZPL width is also a function of the laser excitation 
intensity Plas- 

^ZPl{PlAs) = ^ZPL (0)\/l + Plas/Ps, (5) 


where Ps is the saturation intensity and Tzpl{^) is the 
unsaturated ZPL spectral width. 

From Eqs. m-m it follows that if PpAS ^ Ps 
the temperature is in the range of ultra-low values, then 
the additional broadenings ATe-tunn and ATe-phon can 
be neglected and for each SM it is possible to measure 
the lifetime-limited ZPL spectral width related to Ti. 
Numerous experiments have shown, that these conditions 
are achieved at PpAS < 0.1 — 10.0 W-cm“^ and at T < 
1 — b K (depending on the type of the matrix, see e.g. 
28 . 2^). Thus, SMs can be used as very sensitive probes 
of Ti fluctuations due to the local-field effects. 

If we rule out the sophisticated cases and simulations 
of how structured environments (usually predefined by 
hand) may affect the lifetime of a specific light emitter 
(e.g., recent Sol-liij). the number of correction factors for 
Ti used in practice can be limited to five models 2ol-[24|. 
They all consider the situation when a guest molecule 
emits light to the “continuous medium” of the host ma¬ 
trix, weakly absorbing, isotropic and characterised by n. 
Most generally, the effective Ti is written as [312^: 


Ti(n) = 


To 


n/(n)’ 


(6) 


where tq is the “vacuum” value of the excited state life¬ 
time; n is the refractive index (if precisely, at the wave¬ 
length of the emitted light); function f{n) reflects the 
local-field contributions and includes either squared or 
straight n-dependent coefficient l{n) being the propor¬ 
tionality between the local-field acting on the emitter El 
and the average Maxwellian field Em^ i-e., El = l{n)EM- 


To this date, the choice of l{n) has been among two 
concepts of the local-field. One is the Lorentz local-field 
which is classically calculated under the assumption that 
the field due to polarized molecules inside a small sphere 
centered at the site of a light emitter may be neglected. 
The other is a simplified variant of Onsager’s approach, 
jssi where there is a small empty cavity around the emit¬ 
ter. These cases are often distinguished by the concept of 
interstitial and substitutional guest molecules [25|. The 
most popular models for the lifetime correction are the 
virtual-cavity model 21, 22| based, respectively, on the 
Lorentz local-field: 


/W 


+ 2 


and the real- or empty-cavity model 
Onsager model for the local-field: 


( 7 ) 


based on the 


f{n) = 


3n^ 


2n2 + 1 


(8) 


Another approach is called the fully microscopic model 
0 , in which 


fin) = 


1 + 2 
n 3 


( 9 ) 


Similar to Eq. © it uses the Lorentz field but here the 
l{n) function is not squared. At the same time it misses 
the n factor in the final form of ([6]). The remaining mod¬ 
els imply f{n) = l{n) which brings back the missing n 
and provides agreement with the fully microscopic model 
m in terms that the local-field coefficients l{n) are not 
squared. Thus, for the likely-to-be-interstitial emitters: 


fin) = 



( 10 ) 


while for the likely-to-be-substitutional emitters: 


fin) = 


3n^ 

2n? + l' 


( 11 ) 


Speaking generally, the models in Eqs. © and m are 
based on the field quantization procedures in dielectrics 
and other macroscopic concepts for treating the problem 
of spontaneous emission. The models in Eqs. ©-(HI) ap¬ 
peared from derivation of the Maxwell-Bloch equations 
describing the light-matter interaction using the Heisen¬ 
berg operator formalism for (|9]) and Bogoliubov-Born- 
Green-Kirkwood-Yvon equations for reduced density ma¬ 
trices and correlation operators of material particles and 
modes of the quantized radiation field for ([lOl-idI]). In 
both approaches the field operators were initially written 
for material vacuum providing interactions between the 
guest and host particles. 

Nevertheless, each of these models ©-(HU has found 
“verification” during the analyses of the experimental 
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FIG. 2. The dependence of the excited state lifetime Ti{n) of 
terrylene on the refractive index of its host matrix (squares). 
The curves show the best hts of the Ti(n) dependence us¬ 
ing different models: the virtual-cavity model (red solid); 
the empty-cavity model (brown dash-dot-dot); the fully mi¬ 
croscopic models (green dot); microscopic model in Eq. ()10D 
(black dash-dot) and in Eq. (ITT]) (blue dashed). 

data In all cases the so-called vacuum exci¬ 

tation lifetime tq was subjected to variation during the 
fitting procedure. Depending on the model used for data 
fitting, for each emitting system the preferable Ti{n) de¬ 
pendence was obtained. It must be noted, however, that 
the verifications of the models were based on the data 
from the excited state lifetime measurements for dopants 
in crystals and glasses as well as quantum dots in various 
solutions. It has never been performed for such impurity 
systems as organic molecules in solid matrices with dif¬ 
ferent degree of disorder and, specifically, has never been 
tried on SMs. Here we make the first attempt to test this 
approach using several sets of unique data obtained by 
different research groups. 

In this work we have chosen terrylene (Tr) [38| as the 
object of analysis. It has been widely studied in differ¬ 
ent matrixes by the laser selective spectroscopy meth¬ 
ods, and is one of the most used fiuorophores in SMS. 
For this compound the Ti values were measured for a 
set of solid matrices [s^: polyethylene (PE), polystyrene 
(PS), polyvinyl-hutyral (PVB), polymethyl-methaerylate 
(PMMA), and solid n-hexadeeane (Hex). In the context 
of the present work we also had to find the refractive in¬ 
dices for these matrices. The n values for PE, PS, PVB 
were found in [s^, whereas PMMA and HEX refractive 
indexes were measured by our team using the laboratory 
Abbe refractometer (URL-1, Russia) (see the table in 
Eig. [2]). In Eig. [2] we have plotted the Ti values for Tr 
molecules in all of the above mentioned matrixes against 
the values of n in these systems. The obtained Ti{n) 
dependence was fitted using different models (Eqs. ([7j)- 
dj). The best fit was proved to correspond to the 
virtual-cavity model Eq. ([7]) (red curve in Eig. [2]), with 
the value of tq = 12.1 ns (i.e., Fq = 13.1 MHz). 
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FIG. 3. The distribution of the local values of the refractive 
index in amorphous polyethylene (red circles) and polycrys¬ 
talline naphthalene (blue triangles). The data are normalized 
to the maximum value. The sticks mark the average (for bulk 
media) values of the refractive index of polyethylene (red) [3^ 
and naphthalene (blue) [3^. Black and brown triangles de¬ 
scribe the distributions in naphthalene using the satisfactory 
fit (Eq. (ITOF ) and the fit with Eq. (|8]). The inset shows the 
spectral linewidth distribution for terrylene SMs in polyethy¬ 
lene at 30 mK with Gaussian fit of the data [3^. The same 
for Tr in naphthalene crystal [^ . 

Given the value tq and the best dependence Ti (n) for 
the Tr molecule, we can calculate the local value of n 
at the position of a SM from its lifetime-limited ZPL 
width measured in the matrix involved (Eig. [T]). To do 
this, we have taken the unique experimental data on SM 
ZPL width distributions form Refs. IS El obtained in 
the ETH, Zurich (inset in Eig. [3]). The authors had per¬ 
formed unprecedented and complicated measurements at 
milliKelvin temperatures on the SMS setup equipped 
with the ^He/^He dilution cryostat. At these conditions, 
the broadening contributions AFe-tnnn and ATe-phon 
were negligible. The measurements were performed at 
the laser excitation intensities well below the saturation 
intensity. Thus, in accordance with Eq. (j5j) the obtained 
data were the distributions of lifetime-limited SM ZPL 
widths Fo directly related to SM effective excited state 
lifetimes Ti. These distributions Pro(ro) can be easily 
converted into the distributions of Pti (Pi ) using Eq. ^ 
and then into the distributions of refraction indices P^ (n) 
(Eig. [3]) using the virtual-cavity model (Eq. ([7])) with the 
To = 12.1 ns. 

Considering these distributions of n obtained for the 
polycrystalline and the polymer media one should note 
the following: (a) there are significant fluctuations of the 
refractive index local values in real materials; (b) the 
value of n measured in a bulk sample by classical meth¬ 
ods corresponds to the peak of the P^ (n) obtained in the 
sample. Eig. [3] shows that this is valid for polyethylene 
and naphthalene. This implies the applicability of the 
proposed approach for the probing of the local fluctu¬ 
ations of the refractive index; (c) the local fluctuations 
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of n in the amorphous polymer are substantially greater 
than in the molecular polycrystal. 

If one looks again at Fig. O another satisfactory fit is 
seen (black dash-dot curve), which corresponds to micro¬ 
scopic model (Eq. (p4)|) ). Note that it includes the Lorentz 
local-field factor as in the best fit with Eq. dZD- The other 
three models result in much worse approximations. Yet, 
we have recalculated the distribution of lifetime-limited 
Fq in naphthalene into the distributions of refraction in¬ 
dices Pn(^) using all the rest models. Fig. [3] shows that 
the distributions using the fits for Eqs. m and m are 
subjected to significant broadening and, what is most im¬ 
portant, their peaks are notably shifted from the value 
of n in bulk naphthalene. The other two Pn{n) distribu¬ 
tions not shown are even worse. We treat this as another 
evidence of the correct choice of the model and the va¬ 
lidity of the developed approach. Note that naphthalene 
was absent among the systems (Fig. [2j), in which Ti was 
directly measured. 

To conclude, in our study we proposed a unique ap¬ 
proach for the probing of the local refractive index fluc¬ 
tuations in solids. The method is based on the detection 
of SM ZPLs at conditions allowing lifetime-limited spec¬ 
tral line widths. The great potential of this approach 
is demonstrated. Particularly, simultaneous reconstruc¬ 
tion of SM spatial coordinates with the nanometre ac¬ 
curacy by super-resolution fluorescence microscopy [i^ 
opens the way to perform sub-diffraction refractometry. 
It was found that there are significant fluctuations of the 
local n values in amorphous polymer and molecular poly¬ 
crystalline media. These fluctuations are substantially 
greater in more disordered medium. The peak of the 
distribution Pn{n) corresponds to the value of n, aver¬ 
aged over the bulk sample, which is usually obtained by 
traditional methods. 
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